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Abstract 

We consider the problem of gravitational forces between point par- 
ticles on the branes in a five dimensional (5D) Randall-Sundrum model 
with two branes (at yi and 7/2) and /Z2 symmetry of the fifth di- 
mension. The matter on the branes is viewed as a perturbation on 
the vacuum metric and treated to linear order. In previous work |23j 
it was seen that the trace of the transverse part of the 4D metric on 
the TeV brane, f'^{y2), contributed a Newtonian potential enhanced 
by e^^^^ = 10^^ and thus produced gross disagreement with experi- 
ment. In this work we include a scalar stabilizing field (p and solve 
the coupled Einstein and scalar equations to leading order for the 
case where (p^/M^ is small and the vacuum field 4>o{y) is a decreasing 
function of y. then grows a mass factor e"'^^ where however, fj, 
is suppressed from its natural value, 0{Mpi), by an exponential fac- 
tor e~^^^^''^f^y^ , Xf) > 0. Thus agreement with experiment depends 
on the interplay between the enhancing and decaying exponentials. 
Current data eliminates a significant part of the parameter space, and 
the Randall-Sundrum model will be sensitive to any improvements on 
the tests of the Newtonian force law at smaller distances. An example 
of coupling of the 4> field to the Higgs field is examined and found to 
generally produce very small effects. 



1 Introduction 



Higher dimensional models in particle physics with dimension D > 4 have 
been the subject of much theoretical investigation over the past two decades. 
Higher dimensional theory arises naturally in string/M-theory and is phe- 
nomenologically interesting as they offer the possibility of explaining funda- 
mental features of nature that would not be possible in 4D theory. The 
simplest phenomenology of this type is the 5D Randall- Sundrum model 
(RS1)[T], [2] where the fifth dimension y is compactified with S^/Z2 sym- 
metry so that one can think of space as bounded by two 4D orbifold planes 
(3-branes) at ?/i = and ?/2 = ^rp with boundary conditions at yi and ?/2 to 
enforce the S'^ / Z2 symmetry. With no matter on the branes, the 5D Einstein 
equations have a vacuum solution which preserves 4D Poincaire invariance 
on the branes 

ds^ = e-^^^y\ijdx'dx^ + dy^ (1) 

where 

My) = I3\y\ ; 2/1 - e < 2/ < 2/2 - e ; e > (2) 

and Tjij is the Lorentz metric. Thus if all basic masses are naturally of Planck 
size and the physical world lives on the ?/2 brane, such a structure offers a 
new way of understanding the gauge hierarchy (without undue fine tuning) 
not available in 4D theory. For example, consider a scalar field x the y^ 
brane which we may treat as a perturbation on the vacuum state. The action 
has the form 

S = -\j d'x^ig'^dadjx + m\^) (3) 

where we use the notation /i, z/ = 0,1,2,3,5 and i, j = 0,1,2,3. Letting x' = 
e~^^^X; the theory then takes canonical form with a mass parameter 

rh = e"^^^m (4) 

and the observed mass on the y2 brane would be of TeV size if e~^^^ ~ 10~^^ 
i.e. f3y2 — 35. Thus a Planck size mass travelling on the y2 brane has its 
mass effectively supressed by the strong 5D gravitational forces (much as 
an electron traveling in a solid has its mass modified by the electric fields 
there). The question remains, however, as to whether the 5D theory will 
produce other additional phenomena that would violate known observations 
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on the physical y2 brane. Initial analysis examined whether the Friedmann- 
Robertson- Walker (FRW) cosmology on the y2 brane could be achieved in this 
model [3lll[5l|6l[7l[8l[9l[T0l[IIl[T2l[T3]. This was found to indeed be the case 
provided that in addition to gravity being in the 5D bulk, one must stabilize 
the vacuum metric, which is most easily accomplished by adding a scalar field 
in the bulk, 0(a;*, y) [14J. Then both relativistic and non-relativistic matter 
could be accommodated in the cosmology [71 IS], [T5] . the distance between the 
branes being governed by the density of non-relativistic matter [T3] . 

A second question that has been examined is whether the 5D theory 
correctly reproduces the known gravitational forces between particles. Here 
we treat the matter on the branes as a perturbation to the vacuum metric: 

ds^ = e^'^'^^irjij + hij)dx'dx^ + hi^dydx' + (1 + h5)dy'^ (5) 

There is also a large literature on this subject [l5l[l6l[I7l[T8l[T9l[20l[2ll[22l[23] 
Refs. [15-22] deal only with the gravitational forces on the yi = "Planck" 
brane (where all the masses are of Planck size), and find normal Newtonian 
forces hold between particles on the Planck brane (along with negligibly 
small Kaluza-Klein corrections). In previous work [23], we have examined 
in addition the physically relevant forces on the y2 "TeV" brane, and unlike 
other discussions make sure that the coordinate conditions chosen do not 
lead to bent branes (so that the / Z2 boundary conditions can be correctly 
imposed). To see what occurs for this case, it is convenient to make a 4D 
ADM decomposition of /i^ 



hij = hjj^ + hjj + hi J + hj^i (6) 

where hjj'^ is transverse and traceless {d^hj^'^ = = hl'^'^)and hfj is transverse 
with a non-zero trace f'^{d'^hjj = 0, h]'^ = /^). One may write 

where 

O, ^ # (8) 

(In the above and following, four dimesional indices are raised and lowered 
with the Lorentz metric rjij.) What was found in [23] was that hjj'" gave 
rise to leading order to normal Newtonian forces between particles on the 
Planck or TeV branes. However /-^ gave a Newtonian contribution on the 
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TeV brane that was enhanced by a factor of e^^^^ ~ 10^^, thus producing a 
gross disagreement with experiment. 

None of the analyses discussed above, [l5l|16l[I71[l8l[T9l[20l[2ll[22l[23], 
have included a scalar field (f)Q{x\y) to stabilize the vacuum metric. 
Such a scalar field might produce a mass for the field, thus modifying it 
gravitational potential. In this paper we examine the effects of introducing 
such a stablizing contribution. It is not possible to solve the coupled Einstein 
scalar field equations in closed form, but an iterative solution can be obtained 
when 0q/M| is small (where (poiv) is the vacuum solution and M5 is the 5D 
Planck mass) and (poiv) is a decreasing function. Within this framework we 
find that indeed grows a mass /i but the mass is exponentially suppressed. 
The contribution then appears effectively massless over a distance r < 
which can be anomalously large due to the exponential reduction of /i. 
Thus whether the theory is in agreement with current experimental tests 
of the Newtonian force law at small distances depends on the interplay of 
the amount of suppression of fi compared to the size of the enhancement 
factor e^^^2 of the amplitude of f^, and in fact current experiment strongly 
constrains allowable scalar field models of this type. 

In Sec. 2, we give the choice of coordinate conditions we use and write 
down the Einstein and scalar field equations. In Sees. 3 and 4 we state the 
expansion procedure we use to solve the equations, discuss the solutions of 
the Einstein equations and evaluate at yi and ?/2 (which is what is needed 
to calulate the effect of on particles on the branes). Sec. 5 is devoted 
to the scalar field equations, and in Sec. 6 we calculate the leading effects 
to the fields for a case where a rigorous vacuum solution for 4>o{y) exists. 
In Sec. 7 we calculate the Newtonian forces on particles on the branes. Sec. 
8 examines a (p'^x^ coupling where x is a scalar field (e.g. the Higgs field). 
Conclusions are given in Sec. 9. Appendix A shows that the results obtained 
in Sec. 7 are valid for a general class of models where 0o/M| is small and 
(j)o{y) is a decreasing function of y. Appendix B discusses previous analyses 
of gravitational forces in RSI, and examines why they did not obtain the 
results given here. 
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2 Coordinate Conditions 



The action for our system has the form 

S = J d'xy^l-hllR - A] + y d^x v^[-^5^05^0 - V( 

+ J2 f d'x^/^[Cm^ - V^i^Mv - y^) (9) 

where M5 is the 5D Planck mass, R is the 5D curvature scalar, is the 
Lagrangian for point particles on the Ua brane (a = 1,2), (j){x'^) is the scalar 
field that stabilizes the vacuum metric, V{(f)) is the bulk potential, and Va 
are the brane potentials. We write 

<p{xn = My) + s<Pixn (10) 

where 4>o{y) is the vacuum solution and 6(j) is the perturbation due to matter 
on the branes. The vacuum equations read 

4^" - A" = -^[A + V{<Po)] - ^ E - y'^)U<Po) (11) 

5 5 Q, 

AA- - 4A" = -^[A + V(0o)] - ^0'o^ - ^ E - ^")^"(<^o) (12) 

- 4A'0[, - v'(0o) -J2^^y- ?/")k:(0o) = (13) 

a 

where A' = dA{y)/dy, A" = d'^A/dy^, etc., and V'{(j)o) = dV/dcpQ, etc. 

The bulk and brane potentials are arbitrary except that they must be 
fine tuned to cancel the effects of the bulk cosmological constant A so that 
the net brane cosmological constant vanishes. Most of the analysis can be 
done without specifying V and Va, making use of the field equations Eqs.(ll- 
13). However, to estimate the size of effects, it is useful to have an explicit 
rigorous solution of Eqs. (11-13) and one has been given in [B]. Thus vacuum 
functions 

k2 



My) = 0ie"'^ (15) 
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are the solutions of Eqs. (11-13) for the choice 

Vmy)) = (^^' + 2&/^)<^o-^0o (16) 

Vamy)) = V;(0o(l/a)) + (-l)"260o(l/a)(0o(Z/)-0o(l/a)) 

+ ia\{Hy) - Mvc)? (17) 

where A = —GM^jS'^ (the fine tuning of the cosmological constant), 

VaiMya)) = hir^'mlP - b^liy^)] (18) 

and 7q, are arbitrary constants. We see that the effect of the scalar field is to 
add a term to A(y) of size 0q/M|. Since naturalness implies that all masses 
should be of the same order and comparable to Mpi we will assume 6 ~ 
and b > 0. The gauge hierarchy condition then requires l3y2 ~ 35 [6] so that 
the addition to A{y) is a rapidly decreasing quantity. 

We assume that matter on the branes represent a perturbation to the 
vacuum state and so we solve the full field equations to first order in h^j^^x"^) 
and 50(x"). We begin by reviewing the coordinate conditions we will use in 
the following analysis. The general transformation 

x'^' + e = x'' (19) 

that preserves the /Z2 symmetry with no brane bending is constrained by 

e{x\y{) = Q = e{x\y2) (20) 

As discussed in [23], one may use these to set h^i to zero, but in general it is 
not possible to have h^^ vanish without introducing brane bending. We thus 
assume in the following that 

/i5^ = ; ^ = 0,1, 2, 3 (21) 

There still remains some gauge freedom. Thus under a general transformation 
preserving Eqs. (20, 21), the components of the metric transform to first order 
as ^23j| 

5/155 = 2i% (22) 
6hf = ; 6f = -Q/^e (23) 
5hJ = e'^e (24) 
6{n^h'') = 2e^^a^^^ + 2A'e-2^(e2V),5 (25) 
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where we have decomposed hi and into transverse and longitudinal parts, 
e.g. hi = hj + h\ where d^hj = 0. Eq.(21) requires in addition that 

= e-'^^y^F^x') (26) 
(e^^e^),5 = (27) 

where Ff[x^) is an arbitrary function independent of y. The gauge change 
in 50(a;*) is (to first order) 

5mx')] = <l>'o{y)e{x') (28) 

We see that hfj^ (which contains the Kaluza-Klein modes) is gauge invariant, 
and Eq.(20) implies that f'^ and 6(f) are invariant on the branes at yi and y2- 
Eq.(22) shows explicitly that it is not possible to choose a gauge function 
obeying Eq.(20) that sets /iss to zero everywhere since integrating Eq.(22) 
to try to do this one has only one constant of integration to satisfy the two 
boundary conditions of Eq.(20). As discussed in [23], it is possible however 
to choose a that sets h^^ to zero on each brane 

h55{x\ya) = ; a = 1,2 (29) 

and we will use this gauge in some of the discussions below. 

We conclude this section by recording the field equations and boundary 
conditions. The 5D Einstein equations read 

4^ : ildl - 2A% + \^^U-)h^J = (30) 



g2A 

Ml 



J^T^/ivMy-ya) 



% : ld,v''id,h,, - dih,,) + ^A'd.h, = -^(9,50)0; (31) 



'5 

v'^Rij ■■ ildl - 4A'd,){D^h'^ + f) + e'^D^f + 2A'd,h,, + (32) 

^—n'h, + Ah,{A" - AA") = -^V\ct)M + 

1 Tp^^ 4 2 

^ 5^ 5(y - yo){^ - f'L{<Po)6<P + -h,,VMo)) 

i?55 : (2-95 - A'55)(n'/l^ + f) + ^n'/iss + 2A'95/l55 + (33) 

6 



d'd>Rij 



l\2 



Ml 



3Mi 



V'{<j)o)S<jy 



Ml 



h5{4:A"-AA'^ 
1 Te'^^ 4 2 



+ 



Ml 



3M| 



3M3 



(34) 



In the above T = rj'^-'Tij. The S{y—ya) terms on the right hand side of Eqs.(33) 
and (34) imply that the bulk solutions obey the boundary conditions 



[-l)-+^[d,{n^h'^ + f) + 8A'h,,] 
1 [Te^^- 4^50 + 2/15514] 



3M| 



y=ya 



y=ya 



or equivalently using the vacuum equations and Eq.(40) below 



y=ya 



-,2A 



3M| 



ny) 



y=ya 



(35) 



(36) 



Eqs. (30-34) represent a complete set of Einstein equations. 
The 64> equation reads 

e'^O^ScP - 4A'5cf>' + 5ct>" - y"((/..)50 + h,V'{<Po) + ^d^ia^h"^ + f)<p'. 



-l<p'od5h,, ^J2^(y- ^")[^^55K + v::scp] 



(37) 



with boundary conditions 



14- 



s<i>'{ya) = {-ir^'^[^v:,h,+v:6<j>] 



(38) 



y=ya 
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3 Rj5 Equation 

We consider first the Rjr^ equation. Inserting in Eq.(6) and the orthogonal 
decomposition of hi, Eq.(31) becomes 

ld,{d,f - D'hJ) + ^A'd^h, = ~^d,Scf> (39) 

which can be decomposed into its transverse and longitudinal parts 

d,f + 3A'h,, + 2^S(P = (40) 

dnhj = (41) 

Eq.(41) implies that hj = hj{x'^) is independent of y and one may use the 
remaining gauge freedom of Eqs.(24) and (26) to set hJ to zero, 

hj{x') = (42) 

Each of the terms in Eq.(40) are gauge variant, and it is interesting to see 
how the gauge invariance of the sum arises. Thus using Eqs.(22),(23), and 
(28), the gauge change of the left hand side (Ihs) of Eq.(40) is 

S{lhs) = dy{-6A'e) + GA'^l + (43) 

or 

S{lhs) = {-6A" + 2^-^)e (44) 
Using the vacuum metric equations, Eqs.(ll,12) this reduces to 

s{ihs) ^-^J2^(y- ya)Va{<i>o)e{x\ y) (45) 

5 a 

which vanishes as a consequence of Eq.(20). Thus the gauge invariance of 
Eq.(40) is directly related to the condition that there be no brane bending. 
Eq.(40) allows us to eliminate h^^ in terms of and 54> 
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a relation that holds thoughout the bulk. As mentioned in Sec. 2. it is possible 
to choose a special gauge so that h^^ vanishes on the branes, Eq.(29). In that 
case one has 



y=ya 



y=ya 



(47) 



and one can eliminate d^f^ in terms of 5(f) on the branes. [An alternate 
possibihty is to choose a special gauge such that on the branes 



h. 



55 



y=ya 



(48) 



y=ya 



and then dr^f^ would vanish on the branes. However, in the following we will 
make use of Eq.(29)]. 

4 i?55 and rf^Rij Equations 

Eqs.(32) and (33) give relations between Z"^, , and h^<^. A convenient way 
to analyse these is to first consider the difference Eq.(33)-Eq.(32). Then the 
5{y — Ua) terms cancel and the resulting equation 



3A'a5(n'/i^ + n - e'^n^f + /i55(12A 



/2 



Ml 



(49) 



is valid both in the bulk and on the branes. Eliminating h^^ by Eq.(46) one 
has 



3A' ^ QA'mf''^ 



(50) 



9A'M| A' 



Ml 



One can integrate Eq.(50) to obtain in terms of /"^ and 5(j). The boundary 
conditions Eq.(36) involve precisely the same combination as the l.h.s. of 
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Eq.(50), and since Eq.(50) has been seen to hold on the branes (with no 
^{y ~ Va) singular terms), its solution can be inserted into Eq.(36). In the 
static approximation one has [23] 



T{y^) = -TooiVa) = -e'^^'^-^mj^ir - r{t)) 



(51) 



so that on the boundaries one has 



k' ^2 



20^ 
3 A' Ml' 



'5 



y=ya 



3M| 



(52) 



Finally in the gauge choice of Eq.(47) this reduces to 



2 eT 



□7 



y=ya 



Ml 



-A'Tooiya) - 



(53) 



^[e-^^50(V'(0o)+4A'0[, 



y=yc. 



Note that and 6(f) are gauge invariant on the branes so that Eq.(53) is a 
gauge invariant relation. 

The quantity that governs the Newtonian potential is h^Q, and in the 
static limit this is given on the branes by 



hoo{x\ya) = hl^{x\ya) - -f{x\ya) 



(54) 



Eq.(53) determines the contribution in terms of Tqo and 50. The effect 
of the scalar stabilizing term is to add an additional term, the bracket of 
Eq.(53), and modify the A' factor in the first term, e.g. for the example of 
Eqs.(14) and (15) 



^ 6M| 



(55) 



To examine the effects of these modifications to we consider next the 
field equation. 
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5 5(j) Field Equation 



The 50 field equation, Eq.(37), depends on both /i55 and the combination 
d^iW^h^ + /^). One may ehminate /155 using Eq.(46) and dr:,{U'^h^ + /^) 
by Eq.(50). One gets in this way a rather comphcated equation involving 
only 50 and While is determined on the branes by Eq.(53) (and is 
gauge invariant there), it is gauge variant in the bulk, as is 5<p. However from 
Eqs.(23) and (28), the combination 



1 



tT 



(66) 



is gauge invariant in the bulk. Thus if we eliminate 5(j) in terms of Q, one 
will obtain an equation involving only Q, f^, d^f^, and dlf^ . However, the 
latter three are gauge variant in the bulk, and so gauge invariance implies 
that the coefficients of these three quantities must actually vanish leaving 
an equation involving only the gauge invariant quantity Q. A detailed and 
somewhat lengthy calculation shows that this is indeed the case and the 
equation for Q reduces to the following relatively simple form in the bulk 



e'^n'g + Q" - AA'Q' - V"Q + 



3M| 



Q 

(57) 



Eq.(57) thus gives us an uncoupled equation that determines Q in the bulk. 
One may limit the solution to the branes and impose the boundary conditions 
of Eq.(38). In terms of Q, this reads 



100 



y=ya 



\(-^r\\vLh, + v:Q - \v:(!^)f 



y=ya 



(58) 



One needs 



and using Eqs.(ll-13) 



^ _ '^0 A" 
A' A'-^ 



(59) 



3M| A'^ 



V 



1 (0o)K 
3M| ^'2 



) (60) 
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We interpret the prescription of Eq.(58) to mean the hmit from the bulk as 
y ^ Va and so the last term in Eq.(60) does not contribute to Eq.(58). One 
has then the boundary condition 



V 



li-ir'[[lv:^h,+v:Q-l^v:f] 



(61) 



y=ya 



y=ya 



It is convenient now to make use of the gauge condition Eq.(29) which also 
allows us to eliminate d^f'^ {x\ya) by Eq.(47) (from the Rj^ field equation) 



y=ya 



y=ya 



The Q boundary condition then becomes 



(62) 



1 



V 



A' 3M| A'2 



3M| A 



-1 



KQ--^^v::r 



y=ya 



y=ya 



(63) 



We can also eliminate S(f) in terms of Q in Eq.(53) yielding 

f-11 



y=yc. 



Ml 



-A'TooiVa) + 



-2A 



,1 



U-ir^'<P'oV:-V'-4A'^',){Q 



6 A'-' 



y=ya 



(64) 



6 Leading Order Solutions 

As discussed above the Newtonian potential is obtained from the static ap- 
proximation to hoQ{x\ya) on the branes, given in Eq.(54). hQQ{x\ya) is to 
be obtained by solving Eq.(30) (which is similar to the result of [23] when 
no scalar field was present except for the modification of A{y)). f'^{x\ya) 
on the branes is governed by the coupled equations Eqs.(63) and (64). Since 
Eq.(63) depends on Q', one cannot use it to eliminate Q in Eq.(64) (to ob- 
tain an equation depending only on f^) and one must first solve Eq.(57) for 
Q in the bulk and then insert it into the boundary conditions Eqs.(63) and 
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(64) and use those to determine f'^{x\ya) on the branes. (Thus it is the 
boundary conditions on the branes that couple Q and /^.) The Newtonian 
potential then arises from the 1/r part of hQo(x\ Ua)- Since both Eq.(30) and 
Eq.(57) are decoupled equations, the above analysis is in principle doable. 

An analytic solution of the second order differential equations Eqs.(30) 
and (57) is not possible due to the fact they depend on the complicated 
functions Aiy) and (poiy). We note, however, that the corrections to Aiy) is 
proportional to (Pq/M^. On the 1/2 brane this contains a factor of e~^^^^ ~ 
10"^^ and is very small. If we assume also 0^/M| <C 1, this correction is 
also small on the yi brane and one can consider an iteration scheme based 
on the smallness of (p^/M^. Thus to lowest order, Eq.(30) gives rise to the 
same gravitational potential as in |23] (where no scalar field was present) 

V^'^iVa) = -^ g^'^3 ^[^a^a' + mirris] (65) 

where rria = e'^^^rria are the observed masses on the i/a brane (and rria ~ 
0{Mpi)). Higher order effects are presumably small since they are scaled by 

Eq.(57), which determines Q, shows a similar structure. Thus in the 
example of Eq.(16), V" begins as a constant with (9(0q/M|) corrections and 
the remaining terms have (9(0q/M|), (9(0g/M|), ... corrections. Thus to 
lowest order, in the bulk Q obeys the equation 

e^^^D^g + Q" - 4f3Q' - 7'g = (66) 

where 

7' = r'(0o)]0o=o ; P = [A%o=o (67) 

(and 7^ = 6^ + 46/5 = 0{Mpi) in the model of Eq.(16)). Again, higher 
order corrections should be of (9(0q/M|) and be small. The boundary con- 
ditions Eq.(63) contain corrections to the zeroth order part of O^cp^f^) as 
well as ol{(p'Qf/Ml)f), O{{{(j)'of /Ml)Q). If one neglects all these correc- 
tions, Eq.(66) represents free waves propagating in the bulk (discussed in 
some detail in [6j). They in general give no 1/r Newtonian contribution 
(without extreme fine tuning of parameters). Alternately, one might impose 
a Sommerfeld boundary condition that requires excitations in the bulk to 
arise from matter on the branes. To first order this will occur by including 
the C(0o/^) term, since to zero'th order f'^ is proportional to Too. We thus 
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take as the first order boundary condition 



[Q'-i(-ir'7.<3i 



ya 



where 



7a 



(t>0=0 00 



(68) 



(69) 



Turning to the second boundary equation Eq. (64) one might at first suggest 
that one could ignore the bracket on the r.h.s. as it is (9((0g/M|)/-^), and 
include its effects in by iteration. However, to lowest order, the Too term 
gives a Newtonian piece to ~ 1/r, and if one inserts that into the bracket 
term, one sees that in the static hmit (D^ — > V^) the next approximation 
goes as V~^(l/r) which is infrared divergent. Thus one must include the 
lowest order part of the bracket in the first approximation: 



2 fT 



y=ycx 



Ml 



-/5Too(|/a) + 



(70) 



Ml 



6 /? 



y=ya 



We have kept the Q term on the r.h.s. of Eq.(70) as we will see below it is 
of size 0o/^ as a consequence of Eq.(68). 



7 First Order Solutions for Q and / 

To obtain the part of the Newtonian potential to first order one must 
solve Eq.(66), insert it into the coupled boundary conditions Eqs.(68) and 
(70), and then solve these equations for f^{ya)- To solve Eq.(66) we let 

Qix\y)^ J d'pe^^^Q{p\y) (71) 

and set 

g(p,y) = e2^^i?(p,0 ; e(2/) = ^e^^ (72) 

where = —p^ = {p^Y ~ P^- (Ii^ Eq.(72) m/ (5 is shorthand for (m^//?^)^/^ 
where the branch cut is defined by being real and positive for > 0.) One 
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finds that it! obeys the Bessel equation and so the general solution to Eq.(66) 
is 

Q{p, y) = e'^y[A{p)MO + B{p)MO] (73) 

where u — (A + 7^//3^)^/^, and A and B are constants of integration to be 
determined by Eq.(68). 

It is convenient to introduce the notation 

X, = b/P (74) 

Since by hypothesis b and /3 arc 0{Mpi) we expect A5 = 0{1). We require 
6 > so that (f)o{y) is a decreasing function, and (3 > to achieve the gauge 
hierarchy so that Af, > 0. With this notation 7^//3^ = + 4A{, in the model 
of Eq.(16) so that 

u^2 + Xb>2 (75) 

(Appendix A shows that the above actually represents the leading terms of a 
general model when 0q(7/i)/M| is small and (f)o{y) is a decreasing function.) 
Imposing the boundary conditions Eq.(68) determines A and B to be 

A = ;^[6A^.-i(6)-(^-2-A2)A^.(6)][-^^^^0[,(l/i)r(l/i)^ (76) 



and 



B = -_[e,J._,(6) -{y-2- XMm-'^^^^Uy^)f{y^)\ (77) 

+ -^\iiJ.-Aix) - {y-'i^x^)J.{i^)\{--^{X,-x^)Uv2)f{v2)\ 

where ^1,2 = ^(1/1,2) and D is given by 

D = [6^.-1(6) - (^^ - 2 + Ai) J,(ei)][6iV.-i(6) - (i^ - 2 - A2)iV,(6)] 
-[6iV,_i(a) - (z/ - 2 + Ai)iV,(5)][6J.-i(C2) - (i^ - 2 - A2) J.(6)] (78) 

The boundary condition Eq.(70) reduces to leading order in momentum 
space to 

m^f{yo) = (-l)"^roo(t/.) (79) 



+2^[(-l)°A« - A,]0o[(^J.(O + BN^{i)) + ^e-2^^0or (1/)] 
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where 



Aa ^ g = 0{l) (80) 



To calculate the Newtonian contribution of f^{ya) we take the static limit 
(m^ = — — > — p^) and take the low momentum limit of the right hand side 
(rhs) by limiting — > 0. For ^2 this means we assume 

6 = jef'y^ « 1 (81) 

or we are considering momenta 

p < pe-^y^ fa MpilO"^^ ^ ITeV (82) 

corresponding to distances r > 10~^^cm. (The expansion for ^1 is vahd for 
distances greater than the Planck length) . Thus for all experimental tests of 
Newtonian forces, this expansion is valid. Keeping the leading terms we find 
for f^{y2) that 

AMi2) + BN,{i2) ^ (83) 

JJ bull 
(i/-2 + Ai)A6,, . , , 



\\-\2)(t>if{y2) 



(i/ + 2 + A2)Ab /^ , X WT/ N 

</'l(A6 + Ai)/ (t/i) 

(Z/ + 2 - Ai)A;,^^ X fT/^. \„-2uPy2] 



Gun 

where expanding 1/D gives 



-(A;.-A2)0i/^(2/2)e- 



l/D ^ [1 + ('- + 2-A0(.-2-A2) 2.,.. (84. 

/ (^ + 2 + A2)(z/-2 + Ai)^ ^ (z/-2 + Ai)(z/ + 2 + A2) ' 

Eq.(79) in the static limit then becomes to leading order (the second term of 
Eq. (84) is negligible) 

-ff{y2) = ^^00(1/2) + a^e-(^+^^^)^^M/^(|/2) - f{y,)] (85) 
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where 

_ ?4,.,.(i±M^ (86) 

The e-(2+2A6)/3y2 Eq.(85) comes from the e'^'^y^ of 1/D and the 0o factor 
in Eq.(70) (i.e. 0o = iie-^'^y^). 

The analysis for f'^{yi) is more subtle. Here we find that the numerator 
of AJy{^i) + BN^{^i) terms of Q gives contributions of size e'^^^^ (from e.g. 
N,{i2)Uii) ~ (6/6)") and size e^^^^ (from Ui2)N,{ii)). Multiplying 
by 1/D then gives terms of size e~'^^^y^ and 0(1). The 0{1) term actually 
cancels with the 0{1) term of (f>of^{yi)/6 of Eq.(79), and so one must keep 
the second factor in the bracket of Eq. (84) to get a total result of size e"^''^^^ 
on the rhs: 

-ffiVi) = -J^TM + a'e-^'^''''^^y^{f{y,) - f {y,)) (87) 

Note that Ai does not enter in these leading order results. 

One can now easily solve Eqs.(85) and (87) to get in coordinate space the 
results 

- \f{y.) = i^e--e^^- - l^{m, + m,)(l - e"-) (88) 

and 

1 pT/ \ 1 / , - \ , 1 Gn - -nr 

- o/ yy^l = -Q ("^1 + "^2) + ^ 

6 6 r 6 r 

-i— 6-2^^^^(7711 + m2)(l - e-'^'-) (89) 

O T 

where the Newton constant is given by 
ancOl 

^2 ^ ^2g-(2+2A,)/3,, ^9^) 



^In the "stiff potential" limit, A2 00, Eq.(91) reduces to the mass of the radion 
given in Eq.(6.6) of [26] and Eq.(3.19) of 27J, and for general A2 it is equal to the mass of 
Eq.(30) of [29 . These papers are discussed further in Appendix B. 
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The requirement fi^ > implies A2 > At or A2 < — (4 + A;,). 

In the hmit 0i — > (no scalar field) Eqs.(88) and (89) reduce to the 
results of [23j. The presence of the scalar field does indeed grow a mass for 
the field and so in the limit /xr 3> 1 the remaining 1/r piece in Eq.(88) 
precisely combines with the h'^ff of Eq.(65) to give a total Newton potential 
with Newton constant of Eq.(90) on the TeV brane. (On the Planck brane 
the last term in Eq.(89) is negligible and one gets an additional factor of 
5/3.) However, the factor e^^^^ in the first term of Eq.(88) remains, and the 
mass fi is suppressed by e~^^^^''^^'^^ . One may ask how large r has to get 
so that this anomalous behavior becomes negligible and Newtonian physics 
is reproduced on the TeV brane. As a measure of the effects seen here, we 
assume that the Newtonian force has been measured at the 1% level, so that 
the large dominant term in Eq.(88) implies 

3 



or 



fir > 2py2 - /n(0.03) (93) 

Since we are assuming 0i is small and j3 ~ Mpi we set (pl/M^ = 1/10, 
p = 1.22 X W^GeV and e'^^y^ = 10^2. Eq.(93) then implies 

ji/2g-A6/3y2^ > 4.827 X IQ-^^cm (94) 

where 

Ag(2 + A,)(A2-A,) 
^ = 4 + A. + A2 ^''^ 

Current gravitational force experiments have been done at a separation be- 
tween masses as small as 10/im[30]. Hence we require 

Eq.(96) gives an exclusion contour in the A2 — A;, parameter space. For 
example for A2 = 1, one requires 

\b < 0.67 (97) 

to avoid disagreement with experiment. The general excluded region is 



18 




A2 



Figure 1: Allowed region in the A^ — A2 parameter space for j3 = Mpi. Allowed 
points must lie below the solid lines where /x^ (defined in Eqs.(91) and (86)) 
is correctly positive. The further restriction arising from the measurement 
[5U] of the Newtonian force down to 10/im, Eq.(96), is given by the upper 
dashed curves, the allowed region lying below them. (The lower dashed 
curves represent the restriction if data were improved to verify the Newtonian 
force down to 1/im.) In general Xb and A2 (Eqs.(74)and(80)) are expected to 
be 0(1). 

shown in Fig.l. Thus the absence of any deviation from the Newtonian force 
law already rules out a large amount of parameter space. Fig. 2 plots the mass 
^ as a function of A2 for various values of X^- l//i represents the range of the 
Yukawa-like potential of in Eq.(88) (1/eV = 1.97 x 10" mm). However, 
the range by itself does not govern the strength of the interaction due to 
the anomalous e^^^^ factor in the first term of Eq.(88). Fig. 2 illustrates 
the interplay of these two factors. Thus for A;, small, /x is large and the 
exponential decay suppresses the e^'^^^ enhancement at r > 10/im. However, 
as Xb grows, decreases until finally the e^^^^ factor causes a violation of 
the current data [30]. The lower dashed curve represents the lower limit 
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-10 -5 5 10 

A2 



Figure 2: of Eqs.(86) and (91) as a function of A2 and A5. The lower dashed 
curve is the lower bound for n for the current data r = 10/xm. The upper 
dashed curve would be the bound if experiments were improved by a factor 
of 10. The central region is excluded by the constraint /x^> 0. 

on /J, from Eq.(93) for the current experimental data at r = lO/im. The 

upper dashed curve would be the bound if experiments were improved and 
obtained negative results at r = 1/im. Note from Eq.(86) that /i^ — >■ 00 as 
A2 — >■ — (4 + Aft), and /x^ — > as A2 — >■ A;,, as can be seen in Fig. 2. 

8 Coupling of the Goldberger-Wise Field to 
the Standard Model 

The action of Eq.(9) couples the scalar Goldberger-Wise (GW) field cf) to 
the gravitational field on the branes, which is necessary in order to cancel 
the cosmological constant. This also leads to an indirect non-local coupling 
of to Standard Model matter fields via the coupled boundary conditions 
Eqs.(68) and (70). We consider here the possibility of a direct coupling of 
to Standard Model matter fields. To illustrate what might happen in this 
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case, we consider the simple example of adding to the Lagrangian of Eq.(9) 
the interaction on the TeV brane 



C4>=^9j^f{x\y,)x\x^) (98) 

where x is the scalar field of Eq.(3) (e.g. perhaps the Higgs field) with an 
electroweak sized observed mass fri. Here A < Oil) is a coupling constant 
and we have scaled the interaction by M5 = 0{Mpi) since is a 5D field 
of dimension [mass]^/^ (and in RS models only Planck size masses enter 
fundamental interactions). 

The effect of Eq.(98) on the gravitational fields is simply to add this inter- 
action to the stress tensor. To see the effects on 0, we expand to quadratic 
order in the deviation from the vacuum metric where hij,5(f), and 
treated as first order. Then 

>C<^ = ^^^~"^[& + \v''h,,cf>lx' + 2<f>Mx'] (99) 

where the omitted terms are of cubic or higher order. Wc consider the sit- 
uation for the leading terms of Eqs. (14-18) where (pl/M^ is small. Then 
replacing x by its canonically normalized field x' = e~^^^X) the first term 
becomes an additional mass to the x' field 

>C« = A^m^e-^^/^-x" (100) 

where 

m = M^e-^y^ = e>(TeV) (101) 

We see that the additional mass is suppressed by the factor e^^'^^'^'^ and hence 
is negligible. Similarly, the second term just couples this negligible mass to 
the gravitational field. 

Using Eq.(56) the third term of Eq.(99) becomes 




The last term of Eq.(102) is the same size as the second term in Eq.(99). 
Thus the only term in Eq.(99) that needs further discussion is the first term 
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of Eq.(102) and so we may write 

/ (fP' \ 1/2 p-Ai,/3j/2 

^^-^H^) -^Q^^'^y-)^^"^"^ (103) 

Eq.(103), which is a 4D Lagrangian on the y2 brane imphcs an additional 
term proportional to a 'mass' stress tensor contribution m^x'^ to be added 
to the right hand side of Eq.(68) 




(104) 



In the static limit this term is clearly suppressed compared to the / term in 
Eq.(68) since m^x'^ « C[(Tel^)^] and the factor has the additional e^^^^ 
enhancement. 

It remains therefore to discuss the dynamical effects implied by Eq.(103) 
arising from the Kaluza-Klein (KK) modes of Q. These are a set of discrete 
modes with m„ = (— p^)^^^ obeying the homogeneous equations of Eqs.(66), 
(68), and (70) with no brane matter present (i.e. T^,^ = 0). They presumably 
form a complete set so that in the bulk one can write the general solution of 
Eq.(66) using Eq.(73) as 

Q = e'^y Y}^n{p)UU) + Br,{p)N,{^n)] (105) 

n 

where 

Uy) = y (106) 

The KK masses m„ and the ratio Bn/An are determined by imposing the 
brane boundary conditions. We consider first the boundary conditions. 
For the low m„ modes we may use Eqs.(85) and (87) with —p^ — > and 
Too set to zero: 

mlfiy2) = Afiy2)-fiyi)] (107) 

mlfiyi) = e-'^y^^i'[f{y2)-f{y,)] (108) 
Hence f^ivi) = AV'^(z/2) and 

4-2/32/2 

(m5 - + ^')nm) = ; = J'^, (109) 



22 



Eq.(109) has two solutions. Either 



= /^^ — and f^{y2) = arbitrary (HO) 



or 

,2 / ,.2 a2 „„j 



m^7^/i^-A^ and ^ (2/2) = (HI) 
Since is very small, to leading order Eq.(llO) gives 



m 



/x^ ; f{y,)=e-^Py-f{y,) ■ /^(y^) = arbitrary (112) 



The possibility in Eq.(lll) implies from Eq.(68) 

1 



2 



[g'--(-i)"+VQ] =0 (113) 



The m„ in this case arise from the vanishing of D in Eq.(78) i.e. from the 
algebraic equation 



and Bn/An = - 





aS — 7/3 


= 




a 1(5, where 










-l(Cnl) — 


(^- 


2 + Ai)j,(e„i) 


/3 = UNu 


-1(611) — 


(^- 


-2 + Ai)iV,(ai) 


7 = UJv- 


-l{U) — 


(^- 


2 — \2)Ju{U) 


S = UNu 


-l(Cn2) — 


(^- 


- 2 - X2)N,{U) 



(114) 



(115) 



and 



^na = ^niVa) ] a = 1,2 (116) 

Eq.(112) corresponds to what is commonly called the "radion" mode. How- 
ever, note that this mode arises in Q which is the gauge invariant part of S(j) 
and is not an aspect of the metric components f'^{x\y) or hr,5{x\y) which 
are pure gauge variant in the bulk, as seen in Eqs.(22) and (23). 

In the following we will label the "radion" mode of Eq.(112) by n = 0, 
and the usual KK modes arising from Eq.(114) by n = 1, 2, 3, . . .. We begin 
by examining the KK mode contributions to Eq.(103). 
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A numerical analysis is needed to get accurate solutions of Eq.(114) for 
m„. However, since the lowest m„ is 0{TeV), one has that ^ 1 (provided 
m„ <^ Mpi) while ^n2 > 1, and becomes large rapidly. Thus approximate 
solutions of Eq.(114) can be obtained using the asymptotic forms for the 
Bessel functions: 

tMU - Y + 4^ = (4.^-l) + 2(.-2-A.) ^"^ ^''^^ 

One can qualitatively see where the m„ lie. Aside from exceptional values of 
u = 2 + Xb and A2 (to be discussed below) the coefficient of on the right is 
small, and so r/„ = — z^7r/2 + 7r/4 is close to tt, 27r, Svr, .... As n gets large 
^„ gets large and so the postitions of rin migrates to be close to (2n + l)vr/2. 
The exceptional situation occurs when one fine tunes A2 and Ab so that the 
denominator on the right nearly vanishes, making the postions of rjn start 
near 7r/2, 37r/2 for small n as well. 

In the following analysis, we will assume 

U » 1 (118) 

and calculate the leading terms in the expansion in powers of l/^n2- This 
should be accurate for n > 3, and give a qualitative picture for lower n. We 
also have <^ 1 so that since B^/An — —a/ (5 we can write 

g = J^A^x^F^i/) (119) 

n 

where (^„(y) = rune'^y / (5) 

FM ^ e'^>lMi„) - (|1)^'-^(^±|^);V„(«„)1 (120) 

The second term in Eq.(120) is generally negligible for m„ <^ Mpi unless 
Cniy) is close to ^ni, where it becomes comparable to the first term. 

We turn next to obtain the canonical normalization of the 4D fields 
An{x^)- The Lagrangian generating Eq.(66) is 

J^Q = J d'>x,p^)[-\d,Qgl^d,Q - \{d,Qf - \^'Q'] (121) 

where gij^o) is the vacuum metric. To leading order then 
1 r ry2 

CQ^--jd^xj dye-^fyie^^ydiQrj'^djQ + {d^Qf + (122) 
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so that the kinetic energy is 

Cke = / d^xY^di^nd'Am / e-^^^FMFmiy) (123) 

In general, there appears to be no orthogonahty relation involving quadratic 
integrals of the Fn- However, one may show that the terms with n ^ m are 
1/Cn2 smaller than those with n = m, so that to leading order 

JCke ^-\Y.j d''xdiqnix)d'qnix) (124) 

n 

where 

^n = - ; /x^ = /' dye-''y[FM? (125) 

and g„ are the canonically normalized fields. In general, one must calculate 
the normalization factor numerically (for fixed values of Af, and Ai). However, 
it is possible to get an approximate analytic evaluation by changing variables 
from y to and decomposing the integral into a low part and a high ^„ 
part 

^^l=\f f^e-^^^F^ + 1 f^e-^^yF^ (126) 

We can then estimate the second integral by using the large ^„ asymptotic 
form for the Bcssel functions and the first integral by using the small ^„ form. 
One finds that the first integral is a factor of l/in2 smaller than the second 
so we may write 



2(3 



or 



— o 



2^, P 



/^n2 J JIT TT 

rfen[l + sin(2en-A;.7r)] (128) 



Hence 



/3 . r. 1 



< - Z^Cn2[l - 7^(2 + cos(2en - A.tt) - cos(2 - A^tt))] (129) 



n1 
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Thus to leading order in l/^n2 one has 



7rm„ 

We consider next the mass terms of Eq.(121) 

>C7 = -^f / d'x^/^^Q' (131) 

and 

^'^^-^(^sg)' (132) 

Inserting Eq.(119), and neglecting n ^ m terms gives 

^7 = -h" I r ^y'~''' E (133) 

and hence 

=-Wf y: ,1 r w^^""""^ (134) 

Neglecting the low ^„ part gives 

^--hi''^0j!"'wJ'^^-^ (135) 
which evaluates to 

= -I J E / d'xe-'^y-mnql (136) 
Using Eqs.(67) and (16) this becomes 

>C7^-^(A^ + 4A,) / d^x^^[l + 0(l/U] (137) 



For Eq.(132) we use 



d^Ju = 2(3 J, + i/3^( J.-i - J.+i) (138) 
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Because of the extra factor in the second term of Eq.(138), we keep the 
correction to the leading term in the asymptotic expansion of Ji,±i for 
large Then 



- -jsm{U - ^ - 4)] (139) 



which reduces to 



(140) 



Thus inserting in Eq.(130) gives 

^T = -lf d'xJ2mlqUx)[l + 0{1/U)] (141) 

n 

Thus for large the total mass term for the Qn field is 

C-^jCT + ^^^-I I d'xJ2mlqUx)[l + 0{l/U)] (142) 

n 

and the KK modes of Q have TeV or higher masses. 

We return now finally to Eq.(103) to examine the coupling of the KK 
modes to the SM x' field. Inserting Eqs.(119) and (120) gives to leading 
order 



/,2 , 1/2 ^2 

where m is defined in Eq.(lOl). The expansion of the Bessel function gives 



MU)^sl:^MU-Y + l^ (144) 



and using Eq.(117) this becomes 



MCn2) = + 0{l/en2)] (145) 
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One has then 

- (^) (A) "'fhe-'.^y^ ^ g„(x)x'^ (146) 

Thus since 0?/M| < 1 and /?/M| = C(l), the couphng of the KK modes to 
the SM fields is exponentially suppressed as long as is not anomalously 
smaU, e.g. Ab > 0.1. 

The final contribution to Eq.(103) comes from the "radion" mode 

= e^^y[AUQ + BK{Q] (147) 

where ^o(y) — (nT'o/P)^^^, rrio^ jJi and A and B are given by Eqs.(76)-(78) 
with /^(i/i) and /^(?/2) related now by Eq.(112). Since both and ^2 are 
small (^Q, = {fi/ j3)e^y°') the f'^iyi) parts of Eq.(76) are much smaller than 
the f^{y2) parts. Thus the ratio is {v = 2 + \) 



Similarly, for the B term one gets 



Myi)riy 



U e'(^y-My2)F{y2] 



= e-''(^-^>')f}y^ (149) 



Since as can be seen from Fig. 2, the current tests of Newton's law requires 
Ab < 0.7, again the f^{yi) term is negligible. We thus find 

A = ^{u-2 + Ai)0iiV,(6)Ab(A6 - A2)e2'3j'V^(2/2) (150) 
and using Eq.(84) 



and similarly 

R ^ 

6r(i/) 1/ + 2 + A2 V 2 



S - ^ /"^^""^ (^V^^^/^^-0i/%2) (152) 



28 



Since Co{y) 1 over the entire range of y, we can use the smaU terms of 
Ji, and Ni, in Eq.(147). One finds then in coordinate space that 

Q^{x\ y) = Ce^Pyf{x\ y,) , C ^ VA^^^2(i-a,),.. ^^53) 

Thus f^{x\ 1/2) plays the role of the effective field for the n = mode on the 
4D brane at ^2- 

We next determine the canonical normalization of f^{y2)- Prom Eq.(122) 
the kinetic energy is 

= -1 1 d'xd,f^^^d,f dyC\-'Pye^'^^y (154) 

where we have neglected the non-diagonal terms between the n = and 
higher (KK) modes. Hence writing 

qo{x')^^Jiof{x\y2) (155) 

we have that 

ry2 

t^l^C' dye'^'^'^^^y - ^.^ ^ ^ e^(i+A,)fe (155) 

where qoix^) is the canonically normalized field. Thus 

Qo = [2(1 + h)PY'''e-^^+^''^^y^e''^yqo{x) (157) 
and the contribution of Qo to Eq.(103) is 

4"^ = 2A(^)^/^[2(l + ^^)]'''(]|)'^']i^^o(x)[mV] (158) 

where 

Ms = Mge-^^-^''^^^^ (159) 

Thus the qo coupling is strongly suppressed by the large mass M5 i.e. as A(, 
ranges from 0.1 to its maximum value (from the Newton law data in Fig. 2) 
of \ ^ 0.7, Ms ranges from 5 x lO^GeV to lO^^GeV. 
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We next calculate the mass of the normalized Qo field. From Eq.(122) the 
mass term is 

4°^"^ = -\j d'x j'J dye-'^y[^'Ql + (dM'] (160) 

and performing the y integration gives 

4°)- = I d^x[2['^)iu^ - 2)P^e-'^y^]ql{x) (161) 

Thus the qo mode has an effective mass of /5e-2/3f2 = O(TeV), which combined 
with Eq.(158) implies only a very small effect phenomenologically. 

9 Conclusions 

We have examined here the gravitational forces between point particles on 
the branes in the 5D Randall-Sundrum model with two branes and / Z2 
symmetry (the RSI model). In terms of the orthogonal decomposition of 
the 4D part of the metric of Eq.(5), the static Newtonian forces should arise 
from hoQ = Hqq — /^/3 on the branes y = ya,Oi = 1,2, where Hqq is the 
transverse traceless part of the metric (and also contains the Kaluza-Klein 
corrections) and is the trace of the transverse part of the metric. In order 
to impose the /Z2 boundary conditions correctly, it is necessary that the 
coordinate conditions chosen do not produce brane bending. Thus we assume 
here only that h^i{x\y) = 0,z = 0,1,2,3. While it is not possible to have 
/i55 vanish everywhere, one can assume it vanishes on the branes. hfj^{x'^,y) 
is gauge invariant with respect to the remaining gauge freedom, and f'^ is 
gauge invariant on the branes (and plays the role of the radion). 

Without a scalar field, the amplitude of f^{y2) is enhanced by a factor 
g2/3y2 ^ iq32 jY^ai^ing the theory in serious disagreement with experiment |2 3]. 
In this work we have included a scalar stabilizing field in the bulk 0(x*, y) = 
4>oiy) + ^0) where (poiu) is the vacuum solution and 6(f) responds to matter on 
the branes. The presence of can allow /^(|/2) to grow a mass, suppressing it. 
To examine this possibility we considered the case where 0q/M| was small 
and (poiy) is a decreasing function of y. Then one can obtain analytically 
the leading order corrections. One finds that f^{y2) does indeed grow a 
factor e"'"" but is still enhanced by the e^^^^ factor. Further, the mass fi of 
Eq.(91) is suppressed by the exponential factor e~(^+^'')^^2 where A;, defined in 
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Eq.(74) is positive. Thus whether the RSI model is in agreement with current 
small distance measurements of the Newtonian force law depends upon a 
subtle interplay between the amplitude enhancement and the exponential 
suppression of the mass. Current data eliminates large parts of the parameter 
space. The remaining allowed region is shown in Figs.l and 2. 

The Randall-Sundrum 1 model shows interesting features not intuitively 
expected. Thus the fact that an exponential appears in the metric (a fea- 
ture of the solution of the 5D vacuum Einstein equations) modifies the ideas 
of naturalness. While one would expect that the mass of would scale 
by /3, i.e. ~ /3 = 0{Mpi) (with perhaps a model dependent factor) the 
unexpected feature is the additional (model dependent) exponential factor, 
i.e. IX ~ ^e-(i+-^i))/5y2_ Since exponentials vary rapidly, they radically change 
the 'natural' expectation of the size of fi. Such phenomena are intrinsic to 
the Randall-Sundrum model, since one is using exponentials to create an 
'unnatural' solution of the gauge hierarchy problem. Further, the inverse of 
the very exponentials needed for the gauge hierarchy can enter from metric 
factors appearing in the denominator and do so in the amplitude of i.e. 

~ e^^^2^ ig thus remarkable that the theory can survive the experi- 
mental tests of the Newtonian force law. Improvements of these experiments 
at distances smaller than 10//m will therefore further test the model. Other 
tests of the model which could further reduce the allowed parameter space 
might occur when one introduces couplings of the scalar stabilizing field to 
Standard Model matter on the brane. In Sec. 8 we considered the simple 
example of a coupling proportional to 0^%^ on the TeV brane (where x is 
a scalar field, e.g. the Higgs field). While a full discussion would require a 
numerical analysis, we find from an approximate analytic analysis that the 
effects of such a coupling are very small provided A;, is not anomalously small, 
i-G- ^ 0.1 and obeys the Newton law data constraint < 0.7. 
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A Appendix 



In Sec. 7 we considered a special solution for A and 0o of Eqs. (14-15). We 
show here that this actually represents the leading terms of the more gen- 
eral Eqs. (11-13) when (f)Q{x\yi) is small (i.e. 0o(yi)/M| ^ 1) and 0o is a 
decreasing function of y. 

For the situation considered, we can expand A'{y) and V{(j)o) in a power 
series in 0q 

A'{y)=f3 + 0{<f>l) + ... (A.l) 
V{M = \l'<Pl + 0{<Pl) + . . . (A.2) 

where (5 and 7^ are arbitrary constants. Eq.(13) in the bulk gives then 

0[; -4/30[,- 7^00 + 0(01!) = (A.3) 
To leading order then (since 0o is decreasing) 

00 = 0ie-'^ + . . . ; &>0 (A.4) 

where 

6^ + 4/56 - 7^ = (A.5) 
It is convenient to introduce the parameter 

A. = ^ (A.6) 

As discussed in since 6 > 0, the gauge hierarchy requires /3 > so that 

A, > (A.7) 

and Eq.(A.5) implies 

^ = (A^ + 2)^-4 (A.8) 
Eqs. (11) and (12) imply in the bulk that 

3-4" = % (A.9) 
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and inserting Eq.(A.4) (with the conventional boundary condition A{0) = 0) 
one finds 

showing that the special solutions Eqs.(14) and (15) are the leading terms 
in the more general case. The higher terms in Eqs.(A.4) and (A. 10) are 
determined by the choice of bulk and brane potentials V{(f)) and Vq,(0). 

As discussed in [6], it is still possible to achieve a solution of the gauge 
hierarchy when 6 < and (^o(l/) is an increasing function of y. This situa- 
tion is more complicated than the one treated in this paper since the terms 
{(pKyi) / M^Y Q for example in Eq.(57) might become large and dominate. 
Then an analytic solution as discussed here does not seem possible. 



B Appendix 

We consider here the relation to some of the previous work in this area 
(Refs.[26l[l6lE71E8l[29]). 

In the analysis given here, we have considered gauge transformations that 
produce no brane bending (i.e. maintain the /Z2 boundary conditions). 
Setting /ij5 to zero, the static gravitational force for particles on the branes 
is then governed by 

hoo = hl^-^-f (B.l) 

evaluated on the branes. As seen from Eqs.(20) and (23), hoo{ya) is gauge 
invariant. Thus the formalism describes the static gravitational energy in 
terms of gauge invariant parameters. The usefulness of describing gravita- 
tional energies in gauge invariant language for 4D general relativity is well 
known (e.g. Refs.pH EH]) and we make use of these ideas in the RSI 5D 
model. We first consider Ref . [26] which assumes a metric of the form 

ds'^ = e-^^[7]ij - 2F7]ij]dx'dx^ + (1 + 2G)dy^ (B.2) 

which implies in the notation of Sec. 2 

/i55 = 2G ■ ■f = -6F (B.3) 
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but in addition from the general form of Eq.(6) 




(B.4) 
(B.5) 

(B.6) 



3'' 



Eq.(B.4) implies that the "TT" modes are not being considered. (They do 
indeed contribute to the gravitational forces as seen in Eq.(65).) Eq.(B.5) is 
not a priori obvious, but as seen from Eqs. (41,42), the Rj5 equation combined 
with a gauge choice allows one to set hf to zero. 

The problematical choice is Eq.(B.6). Assuming Eq.(B.6), Ref . |26] argues 
that the didj part of the Rij equation allows one further to set G = 2F i.e. 
/155 = — 2/-^/3. We find in fact from the Rij equation that 



O-'e-'^ildl - 2A'd,]{n-h^ - In + i(/.55 + In = (B.7) 



Thus if Eq.(B.6) holds then indeed /155 = —2/^/3. Thus one must check 
whether Eq.(B.6) can be achieved by a gauge choice. From Eqs. (23) and 
(25) on can easily see that starting in an arbitrary gauge one can in fact 
achieve Eq.(B.6) by the gauge transformation 



where the right hand side is in the arbitrary gauge. However, the Einstein 
equations imply Eq.(36) so that Eq.(20) is violated when Tijja) 7^ 0, i.e. 
matter is on the branes. Eq.(B.6) then implies a gauge with brane bending 
and hence with non-trivial boundary conditions to be imposed. Actually 
Ref. f2B] does the analysis of the radion mass under the assumption that 
there is no matter on the branes, so that the analysis is actually valid. It 
should be noted that the radion (/iss) is not gauge invariant and that while 
(in the stiff potential limit) nir does equal the mass /i of Eq.(91), the former 
refers to the radion mass in the bulk while the latter refers to the short range 
part of the gauge invariant gravitational potential of matter on the branes. 
However with the actual gauge chosen, h^^ = —2/^/3 and the two functions 
are directly related on the branes. 

Ref. [26] does subsequently put Standard Model (SM) matter on the branes 
but still imposes Eq.(B.6) which then is inconsistent with the boundary con- 
ditions. As can be seen from Sec. 2, the coupling of the 5D fields to the 



2e^^n^e = d5{a^h'^--n 



(B.8) 
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SM on the branes actually involves five fields ( hjj'" , h^^, f'^, h^, and S(j)) of 
which one may be eliminated by Eq.(40). The remaining coupled set should 
then be treated quantum mechanically and subject to the matter boundary 
conditions. Thus the quantum field treatment of SM coupling is much more 
complicated than having only a single bulk field F (as in Ref . [26j ) . 

In Ref. [27| the gauge choice of Eqs.(B.5) and (B.6) was also used (referred 
to there as the "Newton gauge" ) but these authors do include matter on the 
branes. They state correctly that this gauge implies brane bending. To 
impose the boundary conditions, Ref.[27j makes a transformation to 'locally' 
Gaussian normal coordinates where h^^ vanishes in the vacinity of the branes. 
This requires a different coordinate frame for each brane. The authors assume 
the local Gaussian coordinates do not have brane bending and impose normal 
boundary conditions there. (A derivation that local Gaussian coordinates are 
in fact free of brane bending was given in |23j . It was also shown there that 
global Gaussian coordinates will generally possess brane bending on one or 
both branes, and so the assumption that this is not the case in Ref. [16] is 
not generally valid.) Thus in principle, the formalism of Ref. |27] should be 
able to deduce the gravitational f orces between particles. They do in fact 
calculate the radion mass (in the stiff potential limit) and show that the 
large r limit recovers the long range Newtonian force. However they do not 
generate the short range corrections to the Newtonian force of Eqs.(88) and 
(89) due to the fact that they make an expansion in powers of D^/m^ (which 
actually sum to the short range corrections) and argue that these correction 
terms are unobservable. They also do not discuss the important exponential 
enhancement of the short range force appearing in Eq.(88) which is what 
makes these short range corrections within striking distance of experiment 
(as seen in Fig.l). 
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